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ABSTRACT 

We discuss the effects of microlensing on the broad emission lines (BELs) of 
QSOs in the light of recent determinations of the size of the broad line region 
(BLR) and its scaling with QSO luminosity. Microlensing by star-sized objects 
can produce significant amplifications in the BEL of some multiple-imaged QSOs, 
and could be very relevant for high-ionization lines. We have identified a group 
of ten gravitational lens systems (~ 30% of the selected sample) in which mi- 
crolensing could be observed. Using standard kinematic models for AGNs, we 
have studied the changes induced in the line profile by a microlens located at 
different positions with respect to the center of the BLR. We found that mi- 
crolensing could produce important effects such as the relative enhancement of 
different parts of the line profile or the displacement of the peak of the line. The 
study of BEL profiles of different ionization in a microlensed QSO image could 
be an alternative method for probing the BLR structure and size. 

Subject headings: cosmology: gravitational lensing — quasars: emission lines 



1. INTRODUCTION 

The change in the continuum flux of quasars by stars or compact objects in intervening 
galaxies (gravitational microlensing) is now a well-established observational phenomenon. 
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Several studies have attempted to resolve the structure of the region generating the op- 
tical and UV continuum by using the microlensing effect as a gravitational telescope (see 
Yonehara 1999, and references therein). However, it has usually been assumed that the 
size of the region generating the broad emission lines (BELs) of quasars is too large (0.1-1 
pc within the framework of the standard model, Rees 1984) to be substantially affected by 
stellar-mass lenses. According to previous studies (Nemiroff 1988; Schneider & Wambsganss 
1990), by comparing one line in different components of a multiple-image quasar, we would 
be able to detect only fractional deviations of the lensed from the unlensed profile. Even 
using optimistic estimates for the microlens masses, these deviations would amount to a 
modest 10-30% range. However, new results concerning the BLR structure and kinematics 
are relevant to these studies, and could substantially change the common view about the 
expected influence of microlensing events on the BEL profiles. A recent study (Wandel, 
Peterson, & Malkan 1999) that compares the BLR size determinations obtained using both 
the reverberation and the photoionization techniques for a sample of 19 low luminosity AGN 
(mainly Seyferts) inferred a size in the range of a few light days to a few light weeks. For 
AGN of greater luminosity (QSOs) Kaspi et al. (2000) derived sizes from the Balmer lines 
in the range from 13 to 514 light days, finding a global scaling of the BLR size as a function 
of the intrinsic luminosity, tblr oc L ar . 

The purpose of this paper is to revisit the influence of microlensing in the BEL by 
incorporating these new results. In §2 we estimate the possibility of microlensing on the 
BLR of different known multiple-image quasars. In §3 we use simple, but not kinematically 
unrealistic, models of BLRs to study qualitatively the effects that microlensing can produce 
on the line profiles. In §4 the observational implications and applications will be discussed. 
Finally, the main conclusions are summarized in §5. 

2. ESTIMATES FOR SOME KNOWN SYSTEMS OF MULTIPLE-IMAGE 

QSOs 

Only extended objects of sizes comparable to or smaller than the Einstein radius asso- 
ciated with the gravitational lens will experience appreciable amplifications (e.g., Schneider, 
Ehlers, & Falco 1992). Thus, in the framework of the AGN standard model where the BLR 
is supposed to have a radius in the 0.1-1 pc range, only massive deflectors could give rise 
to significant amplifications. However, recent results from the MACHO project indicate 
that the most likely microdeflector masses in the Galactic halo are in the range 0.15-0.9 
M Q (Alcock et al. 2000b). In the Galactic bulge the microdeflector masses are in good 
agreement with normal-mass stars (Alcock et al. 2000a). Estimates from the lightcurve of 
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Q 2237+0305 (Wyithe, Webster, & Turner 2000) are also in reasonable agreement with these 
values. Consequently, significant amplifications of the BEL would not be expected accord- 
ing to the standard model. This result has been pointed out by other authors (Nemiroff 
1988; Schneider & Wambsganss 1990), even if they were somewhat optimistic concerning the 
distribution of the microlens mass adopted. 

However, recent research on the BEL seems to indicate that the BLR could be sub- 
stantially smaller than expected in the standard model (Wandel et al. 1999). Based on the 
idea that 'the continuum/emission line cross-correlation function measures the responsivity- 
weighted radius of the BLR (Koratkar & Gaskell, 1991)', Wandel et al. (1999) have obtained 
reliable size measurements. The sizes obtained using this technique (reverberation mapping) 
are consistent with the substantially more extensive but less accurate measurements inferred 
from the photoionization method. The results summarized by Wandel et al. (1999) exhib- 
ited a large scatter (from 1.4 to 107 days) which could be attributed to: i) the different 
size/structure of the BLR in different objects, and ii) the different sizes of the regions asso- 
ciated with emission lines of different degrees of ionization. The latter possibility has been 
neatly exemplified in the case of NGC 5548, in which luminosity- weighted radii ranging from 
~ 2.5 days for the He II A1640 line to the 28 days corresponding to the [C ill] A1909 line 
have been found (Peterson & Wandel 1999). NGC 5548 and the other objects in the sample 
of Wandel et al. (1999) are AGN of relatively low luminosity. However, there is a range 
of luminosities among the lensed QSOs. To take this fact into account we use the results 
from Kaspi et al. (2000), who found a dependence between the BLR size and the intrinsic 
luminosity of the AGN, r blr oc L 7 . 



2.1. Search for candidates 

For a typical lens configuration [z\ = 0.5 and z s — 2), the projected Einstein radius of the 
microlens on the source is r] ~ 2O(M/M )( 1 / 2 ) light-days (for a Q = 0.3 flat cosmology and 
H = 70 km s _1 Mpc -1 ). Objects with r blr ^ Vo are significantly affected by microlensing. 
So the relevant question is whether in the sample of the ~ 70 known gravitational lens systems 
there exist some lensed quasars with BLR radius tblr < Vo- To answer this question we 
can apply the Kaspi relationship tblr oc 

L 0.7 

using; NGC 5548 reference object. We 
adopt two values for the BLR size of NGC 5548 in order to account for its stratification: 
21.2 light-days (Balmer lines, Kaspi et al. 2000) and 2.5 light-days (high ionization, He II 
line, Peterson & Wandel 1999). One can straightforwardly verify that for a z — 2 (1) quasar 
with m v = 24.3 (22.5) would have the same t B lr as NGC 5548, i.e., the microlensing on the 
BLR would be quite pronounced. Multiple-image objects of this apparent magnitude can be 
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detected and, in fact, there are several examples among the currently known gravitational 
lens systems. 

We can refine this rough estimate by taking into account the redshift and apparent 
magnitude of each lensed object. With this aim, in Figure 1 we show contour plots of the 
BLR radius as a function of the source redshift and apparent magnitude using both reference 
values for the BLR of NGC 5548 and the tblr oc L ' 7 law. We have also included in Fig. 1 
the observed redshift-magnitude values corresponding to a sample of 31 QSOs selected from 
the ~ 70 known gravitational lenses with the criteria of having the lens and source redshifts 
and the optical magnitude well determined. We have used the V magnitude of the brightest 
lensed image (in B 1933+507, B 2045+265, MG 0414+0534, PKS 1830-211 and other four 
unlabeled objects we have inferred V from / using V — I — 0.50; in HST 1413+5211 and 
other unlabeled object we have inferred V from R using V — R — 0.11). The data have been 
mainly obtained from the CASTLES web page (http://cfa-www.harvard.edu/castles). For 
the previous and subsequent calculations we have adopted from NED the visible magnitude 
and redshift of NGC 5548, my = 13.3 and z = 0.017. To transform apparent magnitudes 
into intrinsic luminosities we have made use of the equation L = 47i\D 1 2 um S'(l +z) (y(x ^ l \ which 
relates the absolute luminosity of the source, the luminosity distance, D\ um , the apparent 
flux, S, and the spectral index, a, defined by F u oc u~ a . We have considered an Q = 0.3 flat 
cosmology and a value for the spectral index a = 0.5 (e.g., Richards et al. 2001). 

In Fig. 1 we can identify a group of ten systems with magnitude my > 21 that would 
lie in the region R < 100 light-days (Balmer lines) and in the region r < 10 light-days 
(high-ionization lines). These systems are possibly affected by microlensing, with particular 
strength in the case of the high-ionization lines. However, the observed magnitudes of the 
lensed QSOs should be corrected by extinction and gravitational lens amplification. The 
amount of extinction in the gravitational lenses is unknown for most of them, but it can 
be in the range from dust-free lenses to strongly reddened systems (e.g., Falco et al. 1999; 
Muhoz et al. 2001 and references therein). We selected the magnitude from the brightest 
QSO image, which in most of the cases is the less reddened. As a first approximation 
to the expected amount of obscuration we can adopt the mean extinction, (Am) = 0.6, 
derived by Falco, Kochanek, & Munoz (1998) comparing the statistics of radio and optical 
lensed quasars. On the other hand, the quasar source is amplified by the gravitational 
lens and the true luminosity of the unlensed quasar should be calculated by fitting a lens 
model to each system. The exact amplification of each gravitational lens system depends on 
the lens model and on the particular configuration of the system, but an average expected 
amplification between the brightest image and the unlensed source is ~ 1.5 mag (e.g., Lehar 
et al. 2000). Thus, a roughly averaged correction of the combined effects of extinction and 
lens amplification can be made by adding ~ 1 mag to the apparent QSO magnitude. 
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Taking into account the +1 mag shift, the number of possible candidates affected by 
microlensing (at least in the high-ionization lines) is increased to 12 (~ 40% of the sample, 
see Fig. 1). Our selection procedure could include not only intrinsically low-luminosity 
but also reddened objects in the list of candidates. As we have seen, a 1 mag shift due to 
underestimation of the extinction will not substantially modify the set of candidates. This 
moderates the impact of the extinction uncertainties in the selection of candidates, but we 
cannot discard the possibility that our statistics were biased by heavily reddened objects 
of intrinsically high luminosity. This is the case for B 0218+357, which according to the 
rest-frame E(B — V) = 1.52 obtained from Falco et al. (1999) would be an intrinsically high- 
luminosity object reddened by extinction. MG 0414+0534 is also a very reddened object, 
but in this case probably owing to a very red intrinsic spectral distribution (F u oc v~ 9 was 
measured at optical wavelengths by Hewitt et al. 1992). This implies that MG 0414+0534 
could have a high intrinsic rest-frame V luminosity but also indicates that this object is 
probably a radio galaxy rather than a QSO. 

The identification of the candidates available in the literature supports the hypothesis 
that most of them are intrinsically low-luminosity objects, since in most cases the objects 
cannot be clearly classified as bright QSOs but admit alternative classifications as objects 
with lower levels of activity (radio galaxies, underluminous QSOs, starburst galaxies, Seyfert 
2 galaxies, or other types of AGNs). This also implies that we have probably overestimated 
the luminosities of the candidate objects because we have not removed the contribution from 
the galaxy, which is probably far from negligible. 



2.2. Estimate of the amplification in B 1600+434 

It would be very useful to compute the expected amplifications for each object in the 
list of candidates. However, lack of knowledge of the extinction in each system and, to 
some extent, of a precise classification of the source, could affect the results significantly. 
We have selected the gravitational lens discovered by Jackson et al. (1995), B 1600+434, a 
double-imaged typical quasar at z s = 1.59, lensed by an edge-on spiral galaxy at z\ = 0.42, 
in which both the lens amplification and the extinction might be reasonably well known. To 
calculate the amplification induced in B 1600+434 by the lens galaxy we have fitted a singular 
isothermal ellipsoid (SIE) to this double-imaged quasar. Thus, we obtain an amplification 
of 0.96 mag for the brightest quasar image, which has an observed V-band magnitude of 
my = 22.69. For the extinction we have used the result from Falco et al. (1999), Ay = 1.02 
mag. From the intrinsic luminosity obtained after correcting for amplification and extinction 
(my = 22.63) and using the Kaspi et al. relationship, we have estimated radii of 45 and 
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5.3 ligth-days for the Balmer and He II BLRs, respectively (according to the two values of 
reference adopted for NGC 5548). 

To estimate the maximum amplification, // max , we adopt the expression (e.g., Schneider 
et al. 1992): 

^ -, (1) 



' fBLR 



where tblr is the BLR radius, and 770 is given by 



A GM D s D ds 

Vo = a/ 4 —^ (2) 

(D d and D s are the angular diameter distances to the lens and the QSO respectively, and 
D ds is the angular diameter lens-QSO distance). 

If we apply these formulae to B 1600+434, we find that a solar-mass stellar microlens 
would induce a modest, albeit potentially observable, amplification of 1.4 in the Balmer 
lines and a strong amplification of 8.2 in He II. A 0.1 M Q microlens would also induce an 
appreciable amplification of 2.76 in the He II lines. We do not calculate the amplification 
by microlensing on the BEL for the other objects in the list because the exact amount of 
extinction is unknown for most of them, and because we can no longer suppose that the 
F v oc z/ -0 - 5 dependence of the energy distribution is a realistic approach for some of the 
objects. 

In summary, although most of the lensed quasars (> 70% of the total) are intrinsically 
high-luminosity quasars (so that no strong microlensing in the BEL is expected) we found 
that < 30% of the lensed sources are faint enough to be considered as possibly affected by 
microlensing in the BEL. 



3. MICROLENSING EFFECTS ON THE PROFILE OF THE BROAD 

EMISSION LINES 

When an organized velocity field governs the kinematics of the BLR, microlensing can 
give rise to a selective amplification of the emission-line profile (Nemiroff 1988). The shape 
of the line depends on the location of the source with respect to the optical axis (defined 
by the observer and the microlens) and can change with the relative movement between the 
microlens and the BLR. In this section we adopt the kinematic models for the BLR used by 
Robinson (1995) to study the range of profile shapes that exists among AGN in the context 
of a simple parameterization of some basic properties of the BLR such as emissivity and 
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velocity law. Our intention is to study, in this framework, the effects induced in the line 
profile by a microlens in different locations with respect to the center of the BLR. 

We assume that a single microdeflector is affecting the BLR. This supposition adequately 
serves the objectives of the present paper (a more realistic approach based on the existence 
of a random distribution of microdeflectors will be analyzed in a future paper). 

In accordance with Robinson (1995), we consider three different geometries (spherical, 
biconical, and cylindrical) and adopt the following radial dependences for the emissivity and 
magnitude of the velocity: 



= (3) 
and 

We compute the emission-line profile from the expression 

e(r) 5 [A -Ac (l + ^f)] dv , ( 5 ) 



'v 
where 

is the amplification associated with the microlens, r*o ~ (ro, tfo) is the position of the mi- 
crolens, i] is the Einstein radius, and v [t is the projected line-of-sight velocity. 

We adopt inner (r in ) and outer (r blr) radii for the BLR. 

It is convenient to scale the distances to the Einstein radius associated with the mi- 
crolens. In this way, the relevance of the microlens effect can be characterized by the quan- 
tity 

r BLR /_s 
^BLR = • (7) 

We consider two values (1 and 4) for tblr- We assume throughout that r in = 0.1 r B LR- 

To study the effects produced by the relative off-centering between the microlens and 
the BLR, we compute line profiles corresponding to a grid of displacements of the microlens 
relative to the center of the BLR. We consider 25 positions in the positive XY quadrant 
ranging from to f blr in both the X and Y axes (Fig. 2). 

In the Appendix we collect all the formal development and formulae and concentrate 
on the results in the following sections. 



- 8- 



3.1. SPHERICAL SHELL 

In the case of a spherical ensemble of emitters flowing radially, the projected line-of-sight 
velocity corresponding to a emitter at a position (r, 6) is given by 

v,, = v °(i~) cos9 ^ >0 )- ( 8 ) 

We have used the same notation as Robinson (1995) for the parameters related to the 
relative velocity, 

^0 "max 

and we will refer to the line parameter which defines the line shape: 

, = 2±*-l. (10) 
V 

The line profile (see Appendix A.l) is obtained by integrating 

2 / s/3+2-p 

Jo'CkniJ \k) Ux,r,<p)] f=0 drd<f r Iim < r BLR (n) 



F s = < x ov 



r Um > r B LR 



where r Um = r in (i/x m ) 1/p , x m = v /v max = (r- m /r BLR ) p , and [/i(x, r, (p)]f=o is given by Eq. 
6, inserting Eq. A8 into Eq. A10. 

Following Robinson (1995), we take 77 as the parameter defining the emission-line profile. 
For the spherical case we consider two values, a) rj = 2 and b) 77 = —0.25, which in absence of 
microlensing would correspond to concave and convex profile wings, respectively. In absence 
of microlensing, only this parameter is needed to characterize the line profile. In the presence 
of microlensing we need also to fix another parameter. We select the emissivity exponent 
(3 = —1.5 (see Eq. 3). The exponent of the velocity law, p, is then obtained from Equation 
10. Thus, for case a) we have p = 0.5 and for case b) p — 2. 

In Figs. 3, 4 we present the grids of profiles for tblr = 1,4 for case a) and in Figs. 5, 6 
for case b). The influence of microlensing on the line profiles would be observable when the 
microlens is centered on the BLR, and also in many other cases in which the microlens is 
off-centered. The displacement of a microlens across the BLR would induce changes in the 
relative strength of different parts of the line profile, relatively enhancing the wings when the 
microlens is centered on the BLR and the core when it is sited in the outer parts. However, 
no relative enhancements between the blue and red parts of the profile would appear due to 
the high symmetry of the spherical model. 
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As mentioned previously, one notable property of Robinson's models is that the profile 
shape depends only on rj. This implies that from the study of the line profile it is not possible 
to derive direct information on the velocity field or the emissivity law. Could the presence 
of microlensing break this degeneracy for a spherical shell? Such questions as the existence 
of a monotonically increasing or decreasing dependence with radial distance of the velocity 
field may have a formal answer. In principle, it would be possible to invert the line-profile 
expression (Eq. 11) obtained for several different positions of the microlens to recuperate, 
making suitable suppositions, the law for radial velocities. However, it is not easy to decide 
on direct observational criteria to carry out this study. 

To illustrate the difficulties in deriving information from the microlensed line profile in 
the spherical case, notice that not even a simple question like the existence of outflow or 
inflow can be answered. (For the spherical case there is always the same contribution of 
receding and approaching emitters along the line of sight which would undergo the same 
magnification.) 



3.2. BICONICAL SHELL 

Much observational evidence (Zheng, Binette, & Sulentic 1990; Marziani, Calvani, & Su- 
lentic 1992) and theoretical work supports the idea that the flow of emitting gas is anisotropic, 
preferentially confined to a pair of oppositely directed cones. In this model we need three 
polar coordinates (r, 9, ip) measured with respect to the cone axis to express the projected 
line-of-sight velocity corresponding to an emitter, 

v " = v °{i~) z (p>0)j (12a) 

£ = sin # sin sin i + cos # cos i, (12b) 

where i is the inclination of the cones axis with respect to the line of sight. 

This model can give rise to a variety of line profiles (see Figure 5 of Robinson 1995). 
We are going to consider the two limiting cases: i = 0° and i = 90°. We will adopt the cone 
half-angle 9 C = 30°, and (3 = -1.5. 
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3.2.1. i = 0° 

The line profile (see Appendix A. 2.1) is obtained by integrating 

Fx= \ ^ C ImJZZS R] {^Y + P r ' v)]f=odrd V (r lim < r BLR and r sup > r in ) 
[ in the other cases 

(13) 

where r nm = r in (|:r|/:r m ) 1/p , r sup = r B LR (|x|/ cos6> c ) 1/p , and \ji(x,r,<f)] f=0 is given by Eq. 6, 
inserting Eq. A8 into Eq. A10, with x m = v /v max = (r in /r B LR.) p - 

In Figs. 7, 8 we present the grids of profiles for f blr = 1,4 corresponding to 77 = 2, and 
in Figs. 9, 10 the grids corresponding to rj = —0.25. In both cases we obtain two-peaked 
profiles. As in the case of the sphere, the influence of microlensing when the microdeflector 
is centered with the BLR is very noticeable. However, in the biconical case, the change in 
the relative enhancements of different parts of the line profile caused by the displacement of 
the microlens with respect to the center of the BLR is more noticeable than in the spherical 
one. 



3.2.2. i = 90° 



The line profile (see Appendix A. 2. 2) is obtained by integrating 

C Min[7r-6» lim ,7r 



»*BLR 



Fx = 



Max[r lim ,r in 



f(x,r,6) d6 dr, 



where 



f(x,r,0) 



Aofo 



P sing [fi + (x,r,e)+fi-(x,r,e)] f=0 



(9 > 9 hm ) , 







(14) 



(15) 



with 6> lim = arcsin[(|x|/a; m )(r/r in ) p ], x m = v /v B 
given by Eq. 6, inserting Eq. A20 into Eq. A22. 



in other cases, 
(^in/^BLR) p , and [n±{x,r,9)] f=0 is 



In Figs. 11, 12 we present the grids of profiles for tblr = 1,4 corresponding to rj = 2 and 
in Figs. 13, 14 the grids corresponding to rj = —0.25. In this case, i = 90°, we obtain single- 
peaked profiles, convex for rj = 2, and concave for rj = —0.25. The effects of microlensing 
are very strong, even for the case rj = tblr/4 (Figs. 12 and 14). 
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For the limiting geometries considered here (i = 0°, 90°) there are no asymmetries 
induced by the microlensing in the line profile. However, for an arbitrary inclination, an 
off-centered microlens should induce relative enhancements of the blue and red parts. Due 
to the loss of symmetry of the projected velocity field with respect to that of the spherical 
case, the inversion of the profile equation to study the velocity field should be easier. For 
instance, in the case of small cone aperture, the crossing of a microlens along the biconical 
axis would serve to virtually map the radial dependence of the velocity field. 



3.3. CYLINDRICAL SHELL 

Rotation in a plane disk has been often considered in relation to the kinematics of the 
BLR and is typically characterized by the presence of a central dip (Mathews 1982) in the 
line profile arising from the finite extension of the outer radius of the disk. This feature is 
not usually found in the observed line profiles, but there are several ways to avoid it in the 
models (see Robinson 1995). Nevertheless, the existence of rotation is strongly supported 
by recent observational work (Peterson & Wandel 1999). 



3.3.1. Keplerian disk 
The Keplerian disk is a particular case of the cylindrical disk with velocity field, 

( r v 

Vu — v o I — I cosy?sini (p — — 0.5), (16) 



r- r 



where r and ip are polar coordinates of an emitter in the disk. 
The line profile is given by (see Appendix A. 3.1): 



F x 



e r in c rMin[r lim ,r BL R] /_r_ X/3+1 P [n+(x,r)+^(x,r)] f=0 , / x 
Aovosini Jr in \r in J I f " ~- 72 I'lim^'mj 



M^r\ ( i7 ) 

(r lim < r in ) 



where r Um = r in (\x\/ sini) 1 ^, x m = (r in /r BLR ) p and [fj,±(x, r)] /=0 is given by Eq. 6, inserting 
Eqs. A27 and A28 into Eq. A30. 

In Figs. 15, 16 we present the grid of profiles for tblr =1,4 with % = 45°, (3 = —1.5, 
and r] = (P + 2)/p— 1. The most noticeable feature associated with the Keplerian case (and 
with the cylindrical case in general) is the presence of strong asymmetries in the line profiles 
induced by the microlensing. 



- 12 - 



3.3.2. Modified Keplerian disk 

An easy way of generating a single-peaked profile in the cylindrical case is to modify the 
velocity field increasing the contribution of low velocities. For our purpose we will adopt 

v n — v(r) cos</?sini (p = —0.5) (18a) 

and 

v(r) = v ( [_~J™ ) =v u(r), (18b) 

where r and ip are polar coordinates of an emitter in the disk. 
In this case, the line profile is given by (see Appendix A. 3. 2), 

,;;.-.,,< rMin[r lim ,r BL R] ( r \ [n + (x,r)+^(x,r)] f = 



, X v sin i Jr in Uin/ , , L , x n2 V llm m > /1n x 

Fx = \ ^^-{-^f^i) (19) 

(r lim < r in ) 

where r Um = (|x|/ sini)" 1/p (l/r in - 1/V BL r) + I/^blr] , x m = r in /r BL R, and [/i±(x, r)] /=0 
is given by Eq. 6, inserting Eqs. A27 and A28 into Eq. A30. 

In Figs. 17, 18 we present the grid of profiles for f blr =1)4 with i = 45°and (3 = —1.5. 
In addition to the asymmetries, the most interesting effect of microlensing on the line profiles 
corresponding to the Keplerian modified velocity field is the displacement of the line peak 
with respect to the centroid of the non-microlensed line profile, Ax p , (Fig. 19), which, 
independently of the mass considered (r blr = 1,4), can be of the order of as much as 
~ FWZI/4. 



4. DISCUSSION 

If the models and assumptions of the previous sections constitute a good description of 
BLR microlensing, we could conclude that the effects of this phenomenon on the line profiles 
should be not only noticeable but also easily detectable in some lens systems. 

The experimental situation, however, is more complex. In first place, the BELs are 
blended with the narrow emission lines, which come from the much more extended narrow- 
line region (NLR), which would be not affected by microlensing. In fact, it is likely that 
diverse transition regions between the BLR and the NLR could also contribute to the core of 
the line profile. Second, the compactness of the lens systems makes observation very difficult. 
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For instance, in the case (B 1600+434) in which we have obtained a realistic estimate for 
microlensing amplification there are two compact images, the lens galaxy (an edge-on spiral), 
and some additional extended emission, all within a separation of ~ 2". This is a major 
setback when trying to obtain individual spectra, and only modern spectroscopic techniques 
(2D spectroscopy, Mediavilla et al. 1998 and references therein) or observation from space 
avoid the problems induced by source blending and differential atmospheric refraction. In 
any case, to detect microlensing we should obtain spectra (preferably of high-ionization lines) 
with a high S/N ratio in the wings of the line profile, where the contribution from the BLR 
would be dominant. 

The emissivity is another important parameter that can affect the detection of mi- 
crolensing in BELs. In Fig. 20 we present the ratio of the amplified to unamplified line 
profiles corresponding to the modified Keplerian case, where we have changed the emissivity 
parameter, (3. For a highly concentrated BLR ((3 = —4) the effect is stronger when the 
BLR is almost centered on the microlens and weak for larger displacements of the microlens 
with respect to the center of the BLR. But when the emissivity is constant (for a disk of 
uniform brightness, j3 = 0) the effects of microlensing remain noticeable for all displacements 
considered and are hence more likely to be observed. 

In Fig. 1 we have labeled the gravitational lens systems in which significant BEL 
microlensing could be detected (30% of the total). However, in other systems a modest 
effect could be detectable by looking at high-ionization lines. This is important because 
other questions than the amplification, such as the crossing time of the microlens across 
the BLR or the frequency of the events, can lead us to study a gravitational lens system in 
particular. For instance, in a favorable case it would be possible not only to compare the 
line profile corresponding to microlensed and non-microlensed images but also to observe 
in the line profile of an image changes attributable to the microlens crossing. In the most 
favorable case from this perspective, Q 2237+0305, microlensing events are reported each 
year with crossing times of the order of a year or less. However, this is a bright QSO in 
which only very modest amplifications of 30% or less can be expected in the high-ionization 
lines. This estimation of the amplification is, in any case, subject to changes in the BLR-size 
vs QSO-luminosity relationship, and also to the expected intrinsic dispersion from object to 
object. 

In spite of the last comment, the amplifications associated with the high-ionization lines 
could be high in many cases, and only a strong departure from the assumptions made in this 
study (e.g., a severe underestimate of the extinction, very different behavior of the BLR- 
size/QSO-magnitude relationship, or an unexpectedly low-mass population of microlenses) 
could avoid the detection of microlensing in the BELs by comparing the high-ionization line 
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profiles of a microlensed and a non-microlensed image in the most favorable cases. According 
to this, the study of the incidence of microlensing in BELs could become a tool for studying 
the BLR size and stratification, especially when lines of different ionization are observed. 

5. CONCLUSIONS 

In the light of recent discoveries concerning the BLR size and its scaling with AGN 
luminosity, currently accepted values for microlens masses, and a variety of kinematic models 
for the BLR, we have revisited the influence of microlensing in the BEL. We have computed 
grids of line profiles corresponding to different displacements of the microlens with respect 
to the BLR center. Some results are worth summarizing: 

1. - The global amplification of the BEL induced by microlensing events could be relevant. 
We identify a group of ten gravitational lens systems (about 30% of the total sample) for 
which the microlensing effect could be observable, especially in high-ionization lines. In other 
gravitational lenses the microlensing amplifications would be much more modest. 

2. - Even for relatively small microlenses corresponding to high values of the BLR ra- 
dius/Einstein radius quotient, (tblr/ Vo ~ 4), the effects produced by the differential am- 
plification of the line profile (relative enhancement of different parts of the line profile, line 
asymmetries, displacement of the peak of the line, etc.) would be easily detectable except 
for highly symmetric velocity fields. The displacement of the peak of the line profile caused 
by microlensing is especially interesting, since it could otherwise induce inexplicable redshift 
differences between the different images in a gravitational lens system. 

3. - The study of the changes between the BEL profiles corresponding to microlensed 
and non-microlensed images, or among the BEL profiles of lines with different ionization in a 
microlensed image could be useful for probing current ideas about BLR size and stratification. 
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A. MODELING THE BROAD LINE REGION 

Nowadays, there is no generally accepted model for the geometry and dynamics of the 
broad line region. So, we adopt in this paper the discrete cloud model (see Robinson 1995) 
in which it is assumed that the magnitude of the velocity and the volume emissivity may 
be independently specified by a power law. It is also assumed that the clouds or volume 
elements emit line radiation isotropically. 

The amplified line profile produced by the system is obtained from: 

F x = J e(r) 6 [x - A (l + ^)] fi^dV, (Al) 
v 

where the integral is taken over the total volume occupied by the clouds, /i(r) is the relative 
amplification that the microlens produces over the line profile, and v u is the projection of the 
velocity field along the line of the sight. Finally, e(r) is the volume emissivity of the clouds: 

e(r) = e °(£f' (A2) 

r in being the inner radius of the cloud system. 

The geometry of the BLR is unknown, but it seems plausible that its overall structure 
could be characterized by one of three basic symmetries: spherical, conical, or cylindrical. 



A.l. SPHERICAL SHELL 

Let us first consider a spherically symmetric cloud ensemble of outer radius tblr, which 
is characterized by a radial velocity field. The component of the velocity parallel to the line 
of sight is: 

v H = v (—) cos9 (p>0). (A3) 



If the microlens is placed at r ~ (r , y?o) and the Einstein radius is rjo, the relative 
amplification is given by (e.g., Schneider et al. 1992): 

t-* u 2 + 2 If- f | , . . , 

W = -7^=, u = , (A4a) 

|f — fo| = \J (r sin 9 cos ip — ro cos (/? ) 2 + (r sin 9 sin ip — r sin y?o) 2 - (A4b) 
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To compute the integral in Eq. Al, we define: 



/ = A - A 



! + —[ — ] cosfl 



(A5) 



Equation Al, when the integral in the 9 dimension is done and the parameters x 
c (A - A )/(fmaxAo) and x m = v /v max = (r in /r B LR) p are defined, appears as: 



C JSSL,nJ <ry^^^6}^drd V (r lim < r BLR ) 
(r lim > trlr) 



with 



x 



[cos#] /=0 = — — 



and 



[sin 6} f=0 
df" 

d9\ /=0 



[cosfl]^, 



A 



v I r 



[sin^] 



/=o> 



where [//(#, r)]/=o is given by Eq. A4a, inserting Eq. A8 into Eq. A4b as: 



\ r ~ r o|±,/=o = y (r[sin6>] /=0 cosv9 - r cosy? ) 2 + (r[sin 9] f=0 sin (p - r sin^ ) 2 , 
and rii m = r in (\x\/ x m ) 1 / p . The expression of rn m is obtained from the delta condition. 



(A6) 

(A7) 
(A8) 

(A9) 
(A10) 



A. 2. BICONICAL SHELL 



In this case, the broad lines are emitted by a bipolar flow in which clouds moving 
radially outwards are confined to a pair of oppositely directed cones, whose common axis 
passes through the optical axis. In general, an individual cloud measured with respect to 
the cone axis has a projected velocity: 

v 



Vn = V — f (p > 0), 



(Alia) 



£ = sin ^ sin 99 sin i + cos # cos i, (Allb) 

where % represents the inclination of the axis with respect to the line of sight. For this 
system, the amplified line profile is: 



F x = 



e(r) S 



A - Ac 1 + 



ji{f)r sin6d6drd(p, 



(A12) 



where the amplification is given by Eq. A4a, with: 
\f — r | = a/ (rsinOcosip — r cos(/? ) 2 + (r sin 9 sin ip cos i — rcos^sini — r sin(/9 ) 2 - (A13) 
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A.2.1. i = 0° CASE 
The i = 0° case leads to the amplified line profile: 

f2w rMin[r B up,r-BLR.] 



r-Arc rMm[r a up,r B LRj / \ 2 [ft 1 ^] f=0 [ ■ d\ J J / ^ J ^ \ 

Jo JMax[r lim ,r in ] < r ) r fdf ]^ [sm %=o^ (n im < r BL R and r sup > r in ), 







in the other cases, 



(A14) 

where [sin 9]f=o and [df /d9]f =0 are given by Eqs. A8 and A9, \pi(x, r)]f=o is given by Eq. A4a, 
inserting Eq. A8 into Eq. A10, with r\ im = r m (\x\/x m Y^ p , and r sup = r blr i\x\/ cos6 c ) l l p . 
The expressions of r hm and r sup can be inferred from delta conditions. 



A. 2. 2. i = 90° CASE 



In the i = 90° case we consider a projected velocity of 



v,\ = VqI — I sin0sin</? (p > 0). 



(A15) 



Let us now define / as: 

/ = A - A 



1 H — sin sin <£> 

c \ r; •" 



(A16) 



and then, after integrating in the f dimension and adopting that x = c (A — A )/(f m ax^o) 
and x m = v /v max = (r hl /r BLR ) p , Eq. A12 becomes: 



Max[r lim ,r in ] \ Wi 



Min[7r-^ lim ,7r| 



fBLR / Vc Min[7T-W 

I h I 



Tr-e c 



f(x,r,0)d9 dr, 



(A17) 



where 



f(x,r,6) = j [^l-./=o 



j-HFj — — , > PlimJ 

in other cases 



(A18) 



with 6*i im = arcsin[(|:r|/:r m )(r/r in ) p }. The expression of 9\ im is inferred from delta conditions. 
Moreover, 



X ( T \ ^ \ 

[sin99] /=0 = — — — ; — , 
V r in/ sine 



[cosv?]±,/= = ±^1 - [sin^]J =0 , 



(A19) 
(A20) 
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and 



df 

dip 



= -A — (— ^ sin6> [cosy?] ±i/=0 . 



[/i±(x, r)]/=o is given by Eq. A4a, inserting Eq. A20 into Eq. A13 as 

\r- r \±,f=o = \/ (r sin 6» [cosy?] ±J=0 - r cosy? ) 2 + (rcos6> + r siwp ) 2 . 



(A21) 



(A22) 



A.3. CYLINDRICAL SHELL 

The simplest example of cylindrical symmetry is a plane, i.e., a thin disk. We will 
suppose that this disk has uniform thickness, h <^ r- in , and that the angle between its axis 
and the line of sight is i. Finally, any point in the disk is assumed to follow a circular orbit 
about the axis, thus giving a line-of-sight velocity of: 



v u = v(r) cos ip sin i, 

where r and ip are polar coordinates of an emitter in the disk. 
For this system, the amplified line profile is 

F x = J e{r) 5 [a - A (l + ^)] n{f) dS. 



(A23) 



(A24) 



Let us now define / as: 



/ = A - A 



1 + 



v(r) cosy? sin i 



(A25) 



and x = c (A — Ao)/(f m axAo), with v = v max . Then, when the integral in the ip dimension is 
done, Eq. A24 becomes: 



jMHr liM] e(r)r [^(^-(^]/=0 dr (nim > rin) 



Ld^J+./=o 



where 



[cosv?] /=0 = x 



^0 

v(r) sini' 



[sin^± i/=0 = ±sjl - [cos^ =0 , 



and 



dip 



vlr) am 
A [smip\± J=0 . 



(A26) 

(A27) 
(A28) 

(A29) 



±,/=o 
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\p>±(x, r)]f=o is given by Eq. A4a, inserting the expressions of [cos<^]/ =0 and [sin</?]/ =0 (Eqs. 
A27 and A28) into: 

\r- r \±j= = \J (r[cosv?] /=0 cosi - r cos</9 ) 2 + (r[sin y?]±,/=o - r siny9 ) 2 , (A30) 
and rii m is defined by the condition: 

x 

v(r hm ) = v ^-. (A31) 
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Fig. 1. — Contour plots of the BLR radius (in light-days) as a function of source redshift 
and apparent magnitude using both reference values for the BLR of NGC 5548 (R when 
r BLR.(NGC 5548) = 21.2 light-days and r when r BLR (NGC 5548) = 2.5 light-days) and the 
Kaspi et al. relationship r B LR oc L - 7 (for an f2 = 0.3 flat cosmology and H = 70 km s" 1 
Mpc -1 ). Points represent a sample of 31 QSOs whose redshift-magnitude values have been 
observed. No extinction or amplification corrections are taken into account. 
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Fig. 2. — Grid of relative displacement between the microlens and the BLR. The big disc 
represents the BLR. The small discs correspond to the Einstein circle associated with the 
microlens, represented by a point, in the case r] = Tblr/4. For each point (corresponding 
to a displacement of the microlens in the positive quadrant) we compute an emission-line 
profile. 
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Fig. 3. — Spherical model with p = 0.5, (3 = —1.5, and 770 = tblr- On the x-axis we represent 
x = c (A — A ) / (v max Xo) , which varies between —1 and 1. On the y-axis we represent the 
flux. The heavy solid line is the amplified line profile and the solid line is the unamplified 
line profile. 
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Fig. 4. — The same as in Fig. 3, but for r/ = t BL r/4. 
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Fig. 5. — Spherical model with p — 2, /3 — —1.5, and i]q = tblr- On the x-axis we represent 
x = c (A — A ) / (v max Xo) , which varies between —1 and 1. On the y-axis we represent the 
flux. The heavy solid line is the amplified line profile and the solid line is the unamplified 
line profile. 
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Fig. 6. — The same as in Fig. 5, but for r/ = r blr/ 4. 
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Fig. 7. — Biconical model with % = 0°, p = 0.5, (3 = —1.5, and 770 = tblr- On the x-axis 
we represent x = c (A — \q) / (v max \ ) , which varies between —1 and 1. On the y-axis we 
represent the flux. The heavy solid line is the amplified line profile and the solid line is the 
unamplified line profile. 
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Fig. 8. — The same as in Fig. 7, but for r/ = t BL r/4. 
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Fig. 9. — Biconical model with % = 0°, p = 2, (3 = —1.5, and rj = r blr- On the x-axis 
we represent x = c (A — \o) / (v max \ ) , which varies between —1 and 1. On the y-axis we 
represent the flux. The heavy solid line is the amplified line profile and the solid line is the 
unamplified line profile. The figure in the bottom-left corner has been multiplied by a factor 
7.2 for display purposes. 
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Fig. 10. — The same as in Fig. 9, but for 770 = ^blr/4. The figure in the bottom-left corner 
has been multiplied by a factor 6.0 for display purposes. 
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Fig. 11. — Biconical model with % = 90°, p = 0.5, (3 = —1.5, and 770 = tblr- On the x-axis 
we represent x = c (A — \o) / (v max \ ) , which varies between —1 and 1. On the y-axis we 
represent the flux. The heavy solid line is the amplified line profile and the solid line is the 
unamplified line profile. 
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Fig. 12. — The same as in Fig. 11, but for r] = t BL r/4. 
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Fig. 13. — Biconical model with i = 90°, p — 2, (3 — —1.5, and t/q = tblr- On the x-axis 
we represent x = c (A — A ) /(v max A ) , which varies between —1 and 1. On the y-axis we 
represent the flux. The heavy solid line is the amplified line profile and the solid line is the 
unamplified line profile. 
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Fig. 14. — The same as in Fig. 13, but for r] = t BL r/4. 
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Fig. 15. — Model of Keplerian disk with i = 45°, p = —0.5, (3 = —1.5, and r/o = tblr- On 
the x-axis we represent x = c (A — X ) / (v max X ) , which varies between —1 and 1. On the 
y-axis we represent the flux. The heavy solid line is the amplified line profile and the solid 
line is the unamplified line profile. 
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Fig. 16. — The same as in Fig. 15, but for i] = tblr/4. 
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Fig. 17. — Model of modified Keplerian disk with i = 45°, p = —0.5, (3 = —1.5, and 
Vo — r blr- On the x-axis we represent x = c (A — \o) / (v max \ ) , which varies between — 1 
and 1. On the y-axis we represent the flux. The heavy solid line is the amplified line profile 
and the solid line is the unamplified line profile. 
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Fig. 18. — The same as in Fig. 17, but for r) = r blr/ 4. 
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. 19. — Displacements of the line peaks, Ax p , in the modified Keplerian disc model with 
45°, p = —0.5, (3 = —1.5, and r] = r blr/ 4 in the different positions (X , Y Q ). 
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Fig. 20. — Model of modified Keplerian disk with % = 45°, p = —0.5, and rjo = tblr/4. On 
the x-axis we represent x = c (A — A ) / (v max Xo) , which varies between —1 and 1. On the 
y-axis we represent the ratio between the amplified flux and the unamplified flux of the line 
profile. The heavy solid line is for (3 = and the solid line is for (3 = —4. 



